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Abstract 

A new generation of low threshold Imaging Atmospheric Cherenkov 
Telescopes (IACTs) may reach 7-ray energies about 10 GeV with high 
sensitivities and very large collection areas. At these low energies cosmic 
electrons significantly contribute to the telescope background and are in 
principle indistinguishable from 7-rays. In this paper we estimate the elec- 
tron background expected for two configurations of the low energy IACT 
MAGIC. We discuss in particular the reduction of the background caused 
by the geomagnetic field at different locations on the Earth's surface. 
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1 Introduction 



A number of low energy threshold Imaging Atmospheric Cherenkov Telescopes 
are under design or construction and will cover the 7-ray region from 10 GeV 
to 10 TcV. Ground-based instruments [[[] provide high sensitivities and huge 
collections areas of the order of 10 5 m 2 , i.e. orders of magnitude higher than 
those achievable by satellite experiments in the same energy range like GLAST 

In this paper we shall concentrate on the IACT MAGIC |3| || since this 
project aims at the lowest energy threshold. The MAGIC collaboration is cur- 
rently building a 17 meter diameter telescope equipped with a ~ 600 PMT pixel 
camera (MAGIC phase I) which will reach an energy threshold of ~ 30 GeV jl| 
(defined as the peak in the raw differential rate, i.e., before any 7 selection cuts). 
In its second phase the camera will be equipped with high quantum efficiency 
hybrid PMTs (HPD) that will lower the energy threshold down to ~ 10 GeV. 
A third phase may use a camera with Avalanche PhotoDiodes (APDs) of even 
higher efficiency, but we will not consider this for our studies. 

IACTs have to deal with a background of electron- and hadron-initiated at- 
mospheric showers. At very low energies below 50 GeV hadrons rarely generate 
enough Cherenkov light to trigger the telescope. According to [| its contri- 
bution to the background is negligible below 30 GeV. In turn cosmic electrons 
will give rise to showers which are practically indistinguishable from 7-induced 
ones and thus can not be eliminated using traditional 7-hadron discrimination 
methods. They are isotropic and can be suppressed by means of an improved 
angular resolution when the telescope searches for point sources, but they consti- 
tute an irreducible background in extended source searches and diffuse gamma 
flux estimations. 

The purpose of this paper is to estimate the fluxes and detection rates of 
cosmic electrons expected in the new generation IACT MAGIC in its phases I 
and II. In particular we shall carefully consider the effect of the geomagnetic 
field on the electron background observed at ground level. 



2 Geomagnetic Field and Rigidity Cutoffs 

The Geomagnetic Field (GF) is a combination of several magnetic fields pro- 
duced by different sources. More than 90% of the field is generated inside the 
Earth and is known as "main field" or "internal field" . The main field is usually 
described as a dipole (as in the original work of Stormer or more precisely as 
a high order harmonic expansion whose coefficients are based on experimental 
observations and updated every few years, the so-called International Geomag- 
netic Reference Field (IGRF). 

The solar wind interacts with the Earth's field and its effect is described by 
an external field. The currents induced in the magnctosphere and ionosphere 
may create magnetic fields as large of 10% of the main field and are variable 
on a much shorter time scale. Magnetohydrodynamic models are called for in 
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order to describe the whole GF, such as Tsyganenko model [Q or the Integrated 
Space Weather Prediction Model (ISM) || §. The GF bends the cosmic ray 
trajectories preventing low rigidity particles from reaching the Earth's surface. 
(The rigitidy of a particle is defined as pc/ze, where c is the speed of light, p 
is the momentum and ze is the charge of the particle). The minimum allowed 
rigidity is know as rigidity cutoff. In the case of a dipole GF the rigidity cutoffs 
can be calculated analytically In this paper we will draw mainly on an 
expression for the rigidity cutoff (R c ) given by Cookejioj which is based on a 
dipole approximated to the IGRF1980: 

o = 59.4 cos 4 A 

c r 2 [l + (l-cos 3 A sin^ sin</)) 1 /2]2 ^' 

where A is the magnetic latitude, 9 is the zenith angle, <p is the azimuth angle 
measured clockwise from the magnetic north and r is the distance from the 
dipole center (in Earth radius units). Analytical procedures always underesti- 
mate R c due to the effect of the Earth shadow effect. A subset of the trajectories 
intersect the Earth and are physically impossible. This results in a region of 
rigidities some of which are allowed and some forbidden, usually referred to as 
penumbra. The penumbra can be studied by ways of numerical procedures of 
backtracking Jll], |l2). The highest forbidden energy, so-called upper rigidity 



cutoff, is 5-30% higher than R c (l3j |14| . 

R c values in the zenith direction according to equation [l] are tabulated in 
table [3] for some geomagnetic latitudes. MAGIC will actually be built very 
close to a geomagnetic latitude of 30°N. Included in the table are the extreme 
cases of the magnetic pole (90° magnetic latitude), for which all trajectories are 
allowed, and the magnetic equator (0° magnetic latitude) for which the cutoff 
energy reaches a maximum. 



3 The Electron Background 

Electron detection rates 

In the following we will make the assumption that electron-initiated showers 
behave exactly like gamma-initiated ones. 

Figure shows the MAGIC raw collection area (i.e. with no 7 selection 
cuts) at zenith for gamma primaries in the case of the PMT camera (MAGIC 
phase I, figure |l]a) and of the HPD camera (MAGIC phase II, figure |l|b). The 
collection area for 7 primaries was calculated using a Monte Carlo simulation 
of the shower development in the atmosphere and the response of the detector 
to the Cherenkov light produced in the shower (see [Q, |) for details). A trigger 
condition of four next-neighbour pixels with at least 7 photoelectrons in a trigger 
region of 0.8° radius was adopted. Electrons arrive isotropically in the camera, 
so the shower axis was generated parallel to the telescope pointing direction 
("on-axis"), but also with a certain angle 8 with respect to this direction ("off- 
axis") in order to characterize all possible incident directions. In particular 
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(5=0°, 0.4°, 0.8° and 1.2° were simulated. Figure [l] shows the collection area 
(Ss) in MAGIC phase I for these directions. The uncertainties in Ss are always 
below 10%. Figure || shows Ss as a function of S. Ss has been normalized to 
the on-axis collection area (Ss=o) in order to show how the trigger efficiency 
decreases with 5. The dependence of Ss/Sg—o on 6 has been plotted for four 
different energies. For phase II no simulation of off-axis 7 showers was available, 
hence we assume in our calculations that the fractions for phase I hold valid for 
phase II and we can use them to obtain the collection areas off-axis out of the 
collection areas on-axis. 

We can use the raw collection areas along with the electron energy spec- 
trum displayed in figure || (taken from jl5|) in order to estimate the raw rate 
of electron triggers in MAGIC. In doing so we assume that electrons arrive 
isotropically with S <1.5° and that Ss can be calculated from the correlations 
in figure 0. Figure |J shows the electron raw differential rate and integral rate at 
zenith as a function of electron primary energy. The uncertainties in the rates 
are basically determined by the uncertainties in the primary electron spectrum 
which range from ^40% at energies around 10 GeV down to ^35% around 1 
TeV. 

The rate estimated for phase I (2.8±1.2 Hz) is significantly lower than the 
rate quoted in Q (9 Hz). This stems from two reasons. To begin with, we 
have correctly taken into account the decrease of collection area as a function 
of off-axis angle whilst pj assumed that the collection area off-axis could be 
approximated by the on-axis collection area. Besides the authors of B did 
actually not quote the expected electron rate but rather an upper limit to this 
rate based on the upper limit to the electron differential spectrum p6| . The 
electron rates obtained represent a few percent of the hadron rate expected for 
both phases I and II [E El. 

Let us consider now the effect of the geomagnetic rigidity cutoff on the 
electron rates. In table |^ we tabulate the electron fluxes and rates one expectes 
at different locations on Earth as listed in table [|. Also shown is the fraction of 
the polar rate (rate at the magnetic pole). The polar rate is equal to the rate 
that we would obtain neglecting the rigidity cutoff. The table displays the rates 
for both phases of MAGIC pointing to the zenith position. The errors in the 
rates mostly come from the uncertainties in the primary electron flux. It must be 
emphasized that the fractions of the polar rate suffer from no uncertainty since 
the errors in both polar and reduced rates are correlated. For MAGIC phase I 
the effect of R c is small for almost all the locations under consideration. Only at 
the magnetic equator do we observe a 25% reduction in the rate. In contrast for 
MAGIC phase II one expects a noticeable reduction in the rate, especially for 
low magnetic latitudes and as high as a factor of 2 at the geomagnetic equator. 

In all previous calculations we have assumed that there is an abrupt cut 
in the electron energy spectrum at the position of the rigidity cutoff, that is, 
we have ignored the penumbra. We may obtain an upper limit to the effect 
of the penumbra in the electron rate by assuming that all trajectories inside 
this energy region are forbidden and that the upper rigidity cutoff is always 
25% higher than the rigidity cutoff. In this extreme case we obtain further 
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reductions in the electron rate of at most 20%. 
Telescope sensitivity 

We consider now the limits that the electron background pose to the detec- 
tor sensitivity. Contrary to hadron showers, electron showers behave exactly 
like 7 showers. Therefore they cannot be rejected using separation procedures 
based on the shower image characteristics. For point source searches, electrons 
which do not come from the source direction may be rejected within the angular 
resolution capabilities of the detector. But for diffuse or extended sources, no 
rejection is possible. (By "diffuse source" we understand a source which emits 
homogenously inside our camera's field of view, while by "extended source" we 
refer to any non-point source). 

The sensitivity $ of an IACT is generally defined as the 7-ray flux necessary 
to produce a 5<r signal in 50 hours of observation. For simplicity we suppose 
that the 7-ray source has the same spectral index as the electron background. 
The sensitivity for diffuse sources at 5, 10 and 50 GeV for different geomagnetic 
locations can be found in table | A 5% improvement is found when going from 
the magnetic pole to the magnetic equator in MAGIC phase I. The improvement 
in phase II is more considerable. The sensitivity at the equator increases by 13% 
at 10 GeV and by 45% at 5 GeV. 

In the case of point sources we can reject electron showers which do not 
come from the source direction. By using the standard Hillas analysis methods 
(in particular cutting at ALPHA parameter greater than 15°) we can reduce 
the number of electrons by a factor of ~ 6, while keeping ~ 70% of the gammas 
(see page 159 of reference |J. Taking into account the orientation of the shower 
image ("head-tail"), it may be possible to reject 50% of the remaining electrons). 
The polar sensitivity for point source detection in MAGIC phase I is found to 
be <&(> 10 GeV) = 2.4-10~ 10 cm~ 2 s _1 . The proportionality factor between this 
sensitivity and the corresponding diffuse source sensitivity can be applied to 
all the other sensitivities in table || to obtain the point-source sensitivities at 5 
and 10 GeV, and, with good approximation, also at 50 GeV. The effect of the 
geomagnetic field will thus be the same for diffuse and point sources. 

Most of the extended sources constitute intermediate cases between the 
aforementioned point and diffuse sources. In consequence the telescope sen- 
sitivity for a general extended source will be determined by its actual angular 
size, but will always have a value between the sensitivities calculated for diffuse 
and point sources. 

The hadronic background is strongly reduced at energies around 10 GeV 
and the electron sensitivity limits which we have calculated are in the order 
of or above the limit imposed by the hadronic background |Q. Therefore any 
accurate measurement will have to take into account the presence of the electron 
background and the effect of the geomagnetic field. 
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Dependence on zenith angle 

Until now we have restricted our analysis to the zenith position of the IACTs. 
Let us try to predict the geomagnetic cutoff effect at other zenith angles. The 
dependence of the electron rigidity cutoff on zenith angle for a number of az- 
imuth angles and geomagnetic latitudes according to equation (Q) is given by 
figure |{| Whereas the cutoff does not vary strongly with the zenith angle for 
geomagnetic latitudes above 30° (at most a 50% increase from to 90° zenith 
angle in all azimuth directions), it rises by as much as a factor of 3 in the ge- 
omagnetic equator. A given IACT exhibits approximately the same collection 
area and energy threshold for all zenith angles below 30°. Hence let us apply 
the zenith collection areas shown in figure |l| to all zenith angles below 30°. The 
maximal reduction in electron rate takes place in the magnetic equator when 
the telescope points to the east. The maximal enhancement happens also in 
the equator when the telescope points to the west. We have tabulated the rates 
expected for MAGIC phase II in table |^ at geomagnetic latitudes 0° and 30° 
in the west and east directions. The maximal rate reduction is 30% whilst the 
maximal enhancement is 20%. MAGIC phase II sensitivity will thus improve 
slightly when pointing to the west and get slightly worse when pointing to the 
east. The effect of the zenith angle is expected to be smaller in the case of 
MAGIC phase I. 

Conversely the IACT energy threshold increases fast with zenith angle above 
30° (see for example Jl?], |l^, |l9| ) . At a zenith angle of 60° the energy threshold 
is already roughly one order of magnitude higher than the threshold close to 
the zenith. Because the cosmic electron spectrum falls with primary energy 
faster than that of the cosmic ray background or than those of most of the 
predicted 7-ray sources we can expect the electron background to be less and 
less significant as we move to higher zenith angles. In addition the effect of 
the geomagnetic cutoff will decrease with zenith angle since R c does not grow 
as fast as the telescope threshold energy. Electrons triggering the telescope are 
thus never below R c . 

4 Valididy of the Geomagnetic Model 

Equation (|l|) is based on a dipole field. Hence we may consider to improve our 
calculations by applying a more realistic main field such as the IGRF. Never- 
theless because the dipole in equation (|l|) was approximated to the IGRF1980 
value, we do not expect our results to differ considerably from the real values. It 
must be emphasized that no analytical solution for the rigidity cutoff is possible 
by using a field model more complex than a dipole field. 

Besides, we have always neglected in our calculations the effect of the ex- 
ternal geomagnetic field described in section |[ As already mentioned, under 
normal conditions this field contributes up to 10% of the global field intensity. 
The external field, however, is asymmetric due to the interaction of the Earth 
with the solar wind, which gives rise to a long tail extending to the direction 
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almost opposite to that of the Sun (actually in a direction perpendicular to the 
solar wind front, but for simplicity we assume that it extends in a direction 
opposite to the Sun). This introduces an asymmetry in the electron arrival 
direction, since electrons going through the tail are more likely to be deflected 
than electrons coming from the direction of the Sun. In other words, we expect 
the rigidity cutoff to grow in the direction opposite to the Sun. For a given 
arrival direction this is a daily effect. As a consequence we expect the minimum 
electron rate to be reached at midnight in the zenith direction. Since IACTs 
operate only at night, the electron rate will on the whole be further reduced by 
this effect. 

A detailed study would necessitate taking into account the fast variations of 
the external field produced during periods of increased solar activity. In gen- 
eral these periods profoundly complicate the calculation of rigidity cutoffs and 
background electron rates. Careful planning of IACT observations for special 
studies of extended sources or the diffuse gamma background will be necessary 
during these periods. 

5 Conclusions 

We have estimated the electron background expected in both phases of the low- 
energy-threshold IACT MAGIC at different geomagnetic latitudes. The raw 
electron rate expected for MAGIC phase I is approximately 3 Hz in the zenith 
direction at the magnetic pole and drops to 2.5 Hz at the magnetic equator. For 
MAGIC phase II these rates are correspondingly 20 and 10 Hz. The reduction 
with decreasing magnetic latitude reflects the fact that the electron rigidity 
cutoff is minimal at high geomagnetic latitudes and maximal at the equator. 
The sensitivity is also improved when going to lower latitudes by as much as 
50% for $(> 5 GeV) at the equator. The detection rate decreases for some 
azimuth angles reaching a minimum when the telescope points to the west and 
increase for other angles peaking in the east direction. The electron background 
is expected to be very much reduced in high zenith angle observations. 

The conclusion can be drawn that any estimate of the 7-ray flux of an 
extended source or of the isotropic 7-ray flux will have to take into account the 
geomagnetic location of the IACT as well as the specific direction in which the 
instrument is pointing, especially for observations in which the zenith angles are 
less than 30°. 
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Figure 1: Collection areas at zenith for 7 -initiated showers for a) MAGIC phase 
I and b) MAGIC phase II. The solid lines represent the collection area in the 
telescope pointing direction (on-axis showers) while the broken lines correspond 
to different off-axis incident angles. The uncertainties in collection area are 
always below 10%. (From J^J). 
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Figure 3: Cosmic electron differential spectrum (from a number of different 
experiments as summarized in \15y). 
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Figure 4: Differential and integral raw electron rates at the zenith position, as 
expected for both phases of MA GIC. 
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Figure 5: Electron rigidity cutoff as a function of zenith angle for six different 
azimuth angles, namely in north, north-east, east, south-east, west and north- 
west directions. The four curves in each figure correspond to four magnetic 
latitudes: ° (magnetic equator), 30 , 60 ° and 90 ° (magnetic pole). 
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Magnetic latitude 


Rigidity cutoff (GV) 


Magnetic pole 90° 


0. 


75° 


0.1 


60° 


0.9 


45° 


3.7 


30° 


7.3 


15° 


12.9 


Magnetic equator 0° 


14.9 



Table 1: Electron rigidity cutoffs in the zenith direction for several magnetic 
latitudes calculated using equation (Qj. MAGIC will actually be built very close 
to 30 magnetic latitude. 
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Magnetic 


Rate phase I 


Fraction 


Rate phase II 


Fraction 


latitude 


(Hz) 


polar rate 


(Hz) 


polar rate 


Mag. pole 


2.8±1.2 


1.00 


21±8 


1.00 


75° 


2.8±1.2 


1.00 


21±8 


1.00 


60° 


2.7±1.2 


0.96 


20±8 


0.95 


45° 


2.7±1.2 


0.96 


20±8 


0.95 


30° 


2.7±1.2 


0.96 


17±6 


0.80 


15° 


2.5±1.2 


0.89 


11±5 


0.52 


Mag. equator 


2.4±1.0 


0.77 


10±4 


0.48 



Table 2: Electron raw rates and fraction of the polar rate as expected for MAGIC 
phase I and II at different locations of the Earth's surface. The polar rate is 
equal to the rate that we would obtain in the absence of geomagnetic cutoff. 
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Magn. 
lat. 


$(> 5 GeV) 
(cm' 2 s _1 sr _1 ) 


$(> 10 GeV) 
(cm- 2 s -1 sr -1 ) 


$(> 50 GeV) 
(cm- 2 s -1 sr -1 ) 


Phase I 


90° 


4.1-l(r e 


i.o-io- 


4.7-10" 8 


30° 


4.M0- b 


i.o-io- b 


15° 


4.1-10" 


i.o-io- b 


0° 




9.5-10"' 


Phase II 


90° 


1.6-10~ 6 


4.4-10- 7 


3.5-10- 8 


30° 


1.5-10- b 


4.4-10" Y 


15° 


1.2-10- b 


4.4-10-' 


0° 


l.l-10" b 


3.9-10-' 



Table 3: Sensitivities at 5, 10 and 50 GeV for diffuse j-ray sources in both 
phases of MA GIG The table shows the magnetic pole and the magnetic latitudes 
for which some improvement in sensitivity is expected at 5 or 10 GeV. 
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Location 


Azimuth 


0° ZA 


30° ZA 


Lat. 30° 


East 


17±6 


14±5 


West 


18±6 


Lat. 0° 


East 


10±4 


7±3 


West 


12±5 



Tabic 4: Electron raw rates expected in the east and the west directions for 
locations at geomagnetic latitudes ° and 30 °for MAGIC phase II. 
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